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ABSTRACT 


The magnetic field produced by magnetization in Earth's crust and lithosphere can be distinguished from the field produced 
by electric currents in Earth’s core because the spatial magnetic power spectrum of the crustal field differs from that of the 
core field. Theoretical forms for the spectrum of the crustal field are herein derived by treating each magnetic domain in the 
crust as the point source of a dipole field. The geologic null-hypothesis that such moments are uncorrelated is used to obtain 
the magnetic spectrum expected from a randomly magnetized, or unstructured, spherical crust of negligible thickness. This 
simplest spectral form is modified to allow for uniform crustal thickness, eilipsoidality, and the polarization of domains by an 
aperiodically reversing, geocentric axial dipole field from Earth's core. Such spectra are intended to describe the background 
crustal field. Magnetic anomalies due to correlated magnetization within coherent geologic structures may well be superim- 
posed upon this background; yet representing each such anomaly with a single point dipole may lead to similar spectral 
forms. Results from attempts to fit these forms to observational spectra, determined via spherical harmonic analysis of 
MAGSAT data, are summarized in terms of amplitude, source depth, and misfit. 

Each theoretical spectrum reduces to a source factor multiplied by the usual exponential function of spherical harmonic 
degree n due to geometric attenuation with altitude above the source layer. The source factors always vary with n and are 
approximately proportional to if for degrees 12 through 1 20. The theoretical spectra are therefore not directly proportional to 
an exponential function of spherical harmonic degree n. There is no radius at which these spectra are flat, level, or otherwise 
independent of n. 
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1. Introduction 

In the enduring absence of magnetic monopoles, magnetic induction is a solenoidal vector field 
caused by electric currents. Macroscopic magnetic induction, or magnetic flux density, B is often 
partitioned into the magnetic field strength H due to conduction currents and magnetization M due to 
magnetization currents 

B = n 0 (H + M), (la) 

where n 0 represents the scalar vacuum magnetic permeability. This partitioning reflects the difference 
between the extrinsic linear velocity and the intrinsic angular velocity of an elementary electrically charged 
particle. Macroscopic net drift of many such charges gives a conduction current; macroscopic net spin of 
charges gives a magnetization current. 

Here attention is focused upon the slowly varying, quasi-steady magnetic field due to natural sources 
within and near planet Earth. The main source of Earth’s main magnetic field is evidently electric current 
flowing in Earth’s roughly spherical, electrically conducting, ferro-metallic liquid outer core and solid inner 
core. Other sources of the geomagnetic field include weaker currents in the electrically resistive, ferro- 
magnesian silicate and oxide mantle; magnetization in the colder lithosphere and crust; and electric currents 
in the saline hydrosphere, the ionosphere above, and the magnetosphere beyond. In the pervasive presence 
of sources, the geomagnetic field can be represented in terms of poloidal and toroidal vector fields and 
associated scalars (see, e.g.. Backus [1986]); however, the potential field representation remains appropriate 
for those portions of the total field with sources located outside the region of observation. 

1.1 The Geomagnetic Field of Internal Origin 

In an effectively source-free region, the geomagnetic induction Bg is treated as a purely poloidal, 
scaloidal, potential field given by the negative gradient of a scalar potential Vg that satisfies Laplace’s 
equation. This field is mathematically separable into the internal source field Bj nt and the external source 
field B ex j according to the location of the sources relative to a reference surface 

B g = - vv 8 

= - v ( v int + v ext> = B int + B exr (lb) 

If the reference surface is the roughly ellipsoidal base of Earth’s atmosphere, then equation (lb) omits non- 
potential fields from electric currents in the lowest atmosphere (e.g., cloud-to-ground lightening). 

To focus attention upon "the core field" B c and "the crustal field" B x , both externally induced 
internal fields Bj nd (attributed to conduction currents induced within the Earth by time-varying external 
source magnetic fields) and hydrospheric fields (due to conduction currents in the oceans) are omitted 

B int = ' VV int = B c + B x + ( B ind + B h) 


The core field is mainly produced by conduction current in Earth’s core (J c = VxH c ). The crustal field is 
mainly produced by magnetization current (C x = VxM x ) in the portion of the lithosphere lying 
gravitationally above, but thermally below, the Curie isotherm. Within this magnetic crust, part of the 
magnetization is induced by the present core field and is parallel to B c (M ind = xH c for susceptibility x); 
however, thermoremanent, viscous, chemical, and detrital magnetizations also contribute to M x . 
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At time t and position r in geocentric spherical polar coordinates (radius r, colatitude 0, and east 
longitude <|>), the spherical harmonic expansion of zero mean V im is well-known to be 

v { \ v a n+1 n 

V int (r,t) ‘ a 1 H s tg n m (t)cosm<J) + h n m (t)sinm<t>]P m (cos0) , Hdl 

n=l m=0 * 

where P n m is the Schmidt-normalized associated Legendre function of degree n and order m and [g m (t), 
h n (t)] are Gauss coefficients at the reference sphere of radius a (about 6371 km). Gauss coefficients can 
be estimated by a weighted least squares fit to geomagnetic data; this usually requires truncating the sum 
over n at finite degree N (see, e.g., Langel [1987]). 

By (lc), each internal field coefficient is the sum of core and crustal contributions 

Kn m 0) = g n ?<0 + 8n?(t) 

h n m W = h nc(‘> + h n ?(0. (le) 

Mathematical analysis of geomagnetic measurements made on and above Earth’s surface cannot uniquely 
separate core from cntstal field coefficients. Geophysical analysis and interpretation of such measurements 
nonetheless distinguishes core from crustal fields. 

1.2 The Spatial Magnetic Power Spectrum 

Lowes [1966, 1974] and others (see Cain et al. [1989b]) show the mean square magnetic induction 
averaged over a sphere of radius r enclosing its sources to be 

1 2n n oo • 

<[B(r,t)] 2 > - / [B(r,t)] 2 sin0d0d«j> = I R^t), (2a) 


where R n (r,t) is the mean square field configured in harmonics of degree n averaged over the sphere 
3 2n+4 n ? 9 

R n (r,0 = (n + I) (— ) X [g n m (t)] + [h m (t)] . (2h) 

r m=0 

The values of R n (r,t) form the spatial power spectrum of the internal magnetic field. Indeed, the magnetic 
energy density per harmonic degree integrated over the sphere, 2jtr 2 p () -lR n (r,t), is the discrete (line) 
spectrum with dimensions of energy per unit length, hence spatial power or force (SI units of Joules/meter) 
Although the SI (MKS A) units of R n (r,t) are (Tesla) 2 , the R n (r,t) are nonetheless collectively refetred to as 
the spatial magnetic power spectrum. An individual value of R n (r,t) implies the spatial power at degree n 
and is thus called a multipole power: R j is dipole power, R 2 is quadruple power, etc. 

The exponential attenuation of R n (r,t) with n ensures convergence of the sum in (2a) at all r > a if the 
reference spectrum at time t, 

R n = (n + l) X (g n m ) 2 + (h^) 2 = Um R n (r), f2c> 


is of finite power law form (R n < Kn a for finite a). Divergence of a summed reference spectrum X_ R 
merely indicates sources at r > a, such as those due to Earth’s ellipticity and surface topographyTe.g." 
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magnetization of Mount Kilimanjaro). Curiously, crustal sources above the reference sphere imply that 
(Id) alone is not enough to represent the field on Earth’s aspherical surface due to sources within the Earth. 
This particular complication in representing the field does not arise at satellite altitudes. 

1.3 Core and Crustal Spectra 

As core and crustal magnetic fields change with time, numerical values of core and crustal 
coefficients realized during a very long interval of geologic time define two populations of coefficients 
(two ensembles). Statistical properties of these populations can be estimated by analysis of geo- 
paleomagnetic observations and used to test theoretical predictions derived from geophysical hypotheses. 
Such predictions may involve probability distribution functions (PDFs) for the coefficients or functions 
thereof, notably expectation values. An expectation value of a quantity q, denoted { q } , is the integral of q 
multiplied by a PDF for q over all possible q. The correlation between two quantites p and q is the 
normalized covariance f(p- {p})(q- {q})}({p^}{q^))‘^. 

To the extent that core and crustal coefficients are uncorrelated and of zero expectation value, 

i(g n m ) 2 ) = i(g n £ + g n ?) 2 i = ug n *?) 2 } + (2g n ^g n ^) + {(g n T> 2 > 

= Kgn 1 ?) 2 ! + {(SnT> 2 > ' (3a) 

A similar relation holds for {(h n m ) 2 (. To the extent that (3a) holds, the expected total spectrum (R n ) is 
the sum of expected core and crustal spectra, 

{ R n l - {IW + < R nx> ■ (3b) 

Equation (3a) implies (3b), but equation (3b) does not imply (3a) because of the sum over harmonic order 
implicit in (3b). Indeed, (3b) can be an excellent approximation, even if some core and crustal coefficients 
are strongly correlated, provided the 2n+l cross correlations, which may be either positive or negative, tend 
to cancel when summed over order m (e.g., (g n £ « .{ g „7 g( 7} * 0 ). So (3b) is less restrictive than 
(3a). It may thus seem reasonable to suppose that 

R n ** ^nc + R nx - (3c) 

Given (3b) alone, a particular set of cross terms (like 2g n ™ g n *J) for a particular degree n amounts to 2n+l 
samples of a population with zero mean. If these samples are random, then the relative error in (3c) ( I R n - 
R nc ‘ R nx I ^n^ * s ex P e cted to fall off in proportion to (2n+l)'^ 2 by central limit theorem. For 
example, the expected relative error in (3c) would be <20% for n £ 12. 

Expectation operators are linear operators; therefore, in regions where (2a) and (2b) hold, (3b) implies 

{<B 2 >} = {<B c 2 > + <2B C *B X > + <B x 2 >} 

- (<B c 2 >) + {<B x 2 >> . (3d) 

Equation (3d) also holds to the extent that the core and crustal fields are uncorrelated in the mean ; it is less 
restrictive than (3b) due to the implicit sum over degree and, in turn, is much less restrictive than (3a). Of 
course, 2B C *B X varies from place to place and typically dominates the scalar anomaly I B I - I B c I mapped 
by many magnetic surveys because I B c I » I B x I . Provided B c is dominated by low degree, broad scale 
structure and B x is dominated by high degree, narrow scale structure, orthogonality of the spherical 
harmonics ensures that l<2B c «B x >| « < 1 2B C *B X I > on the reference sphere. 
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1.4 The Core Spectrum 

According to the famous interpretation of the spatial magnetic power spectrum derived from 
MAGSAT by Langel & Estes [1982], the core field dominates R n for degrees n < 12 and the crustal field 
dominates R n for n > 16. For degrees 1 < n < 12, and apparently 1 < n < 12, the spatial magnetic power 
spectrum is well described by McLeod’s Rule for the core field [McLeod, 1985; 1994, 1996 equation (20)]. 

In a simplified derivation of McLeod’s Rule for finite degrees n [Voorhies & Conrad, 1996], random 
geomagnetic secular variation (OjB) is induced by fluid motion near the top of a roughly spherical, high 
conductivity core of mean radius c. The spatial power spectrum for such secular variation is expected to be 
that of local differential dipole moments of random orientation and random moment 


a 2n+4 

F nc (c > = (n + IX—) 


m=0 


“ ^ F nc^ c ^ — C n (n + 1/2) (n + 1) , 

(c.f. section 2.2). The variation at the top of the free-stream must, however, diffuse through a thin viscous 
boundary layer separating the free-stream from the electrically resistive mantle (see Appendix A). The 
diffusive time-scales depend upon the effective scale heights (irk^)' 1 for poloidal field harmonics 
within the core, as well as core conductivity o, and are 


*n s ({Rnc>/( F nc » 1/2 = ^(jU ^)* 2 


(see (A7)). The effective scale heights are taken to be directly proportional to the horizontal length scales, 
^ “ k hn* s c fa( n + 1)]' 1/2 , and the consequences for x n 2 verified empirically [Voorhies & Conrad, 

1996]. Combining this relation with (4c) and (4b) yields 

R nc «) - IRnc «c„ - K- ^ («, 


- K’ / (n + 1/2) , (4e) 

which, upon upward continuation, is McLeod’s Rule. For degrees n < 12, the implied contribution to the 
reference spectrum is 

^nc “ « K (n + 1/2) * (c/a) 

Geophysical parameters c = 3.5 Mm and K’ = 5.5xl0 10 nT 2 in (40 are determined empirically via analysis 
of geomagnetic observations. 

Voorhies & Conrad [1996] used McLeod’s Rule and ancillary hypotheses to accurately predict (i) the 
mean radius of Earth’s core; (ii) mean paleomagnetic field intensity; and (iii) the mean rates and mean 
durations of major geomagnetic dipole power excursions and durable reversals of Earth’s axial dipole. The 
predicted core radius is within 1% of the seismologic value (c s - 3480 km). The predicted mean 
paleomagnetic intensity is well within the broad range of experimental mean values. The predicted mean 
reversal frequency (2.26/Ma) is within 3% of the paleomagnetic mean value for the past 84 million years. 
Again, the core field spectrum is well described by McLeod’s Rule. 
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I. 5 A Level Crustal Spectrum? 

By analogy with Lowes [1974], when the R n (r) for some range of n are approximated by the 
exponential R n (r) = A (b*/r) n+ , linear regression through observational values of ln[R n (a)] calculated 
from a geomagnetic field model gives 

ln[R*(a)] = n[In(b*/a) 2 ] + [ln(A*) + 41n(b*/a)] . (5) 

The^ slope of this line implies the radius b* at which the model spectrum R*(r) becomes independent of n 
(or levels off ). If R* is extrapolated to arbitrarily high degree, then the sum over n of R*(r) diverges for r 
< b*; therefore, such extrapolation is not valid below b*. This might be due to failure of the potential field 
representation, so b* might be the minimum radius of a sphere containing the sources. 

The Langel and Estes [1982] fitted two spectra of form (5) to the R n spectrum from their degree 23 
model MGST 10/81 of Magsat data, one for the core field and one for the crustal field. Their core spectrum 
levels off 174 km below c. Their crustal spectrum (R^a) = 37.1(0.974) n nT 2 ) levels off 83 km below a. 

Voorhies [1984] fitted the mean square radial field component alone (Bjj, 2 = (n+l)R /(2n+l) => 

R n /2) from MGST 10/81 . His core spectrum levels off 200 km below c; his crustal spectmnUB* 2 (n) = 

I I. 91(0.9969)" nT ) levels off 9.9 km below a. This may suggest upper crustal, rather than deep 
lithospheric, sources. 

Cain et al. [1989b] used their degree 63 model M07AV6 of Magsat data [Cain et al„ 1989a) to obtain 
a core spectrum that levels off about 76 km below c; a crustal spectrum (R* x (a) = 19.1(0.996)" nT 2 ) that 
levels off 14 km below a; and an estimated noise level of 0.091 nT 2 

The results of these analyses are clearly consistent with lithospheric sources for the internal field of 
degrees 16 or more. More sophisticated crustal magnetic spectra advanced by McLeod and Coleman 

[1980], Jackson [1990, 1994], and McLeod [1994, 1996], however, indicate substantial modification of the 
plain exponential form. 


2. Theoretical Magnetic Spectra for Thin Spherical Shells 

Consider the crustal magnetic field B x = -W x caused by all K magnetic domains within Earth’s 
magnetic crust. These domains may be considered microscopic rock magnetic domains within the mineral 
grains comprising Earth’s crust; however, the formalism is the same if each domain represents a 
macroscopic geologic structure. At altitudes above Earth’s surface exceeding several times the domain 
size, each domain k = 1, 2, 3, ... K produces a field that can be fairly well described by that of a point 
magnetic dipole of moment m k = M k /p 0 at domain position r k . The total crustal potential at observation 
position r is thus treated as the superposition of K offset dipole potentials V x k , 


V x (r) 


K 

Z V x k (r) 
k=l 


K M k 1 

- Z .V, 

k=! 4jt ''•*1*1 


(6a) 


K M k . ! 

X • V 

k _ 1 4n I r - r k I ( 6t> ) 


In (6a) and henceforth "M” represents normalized dipole moment, not magnetization. The SI (MKSA) 
units of M are Tesla meter 3 (Tm 3 ). In (6b), V* represents the gradient operator written in terms of, and 

acting on, r k coordinates (r k , 6 k , 4^). The spherical harmonic expansions for the V k are used to rewrite 
(6a,b) above the crust as 
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a n+1 


V?c(r) = k^i a n ^i ( ~ } m f 0 tA " kCOSm4 * + B n k sinm WP n m (cose) . (( 

To find the coefficients (A n JP, B n ^>) for moment M k , recall that 
1 « n (r k ) n 

(TVq = n =0 m=0 7^1 [C n m(0 ^ )C n m < 0 k%) + S n ”C©4 »S n m (B k ^)] , C 

^emL" noT?r (COS0) and Snm S Sinm<t,P n m ( cos0 ) <«*. e*~ ^kson [1975, equation (3.70)] , 
gel [1987, (195)]. Using (7) to evaluate the gradient in (6b), and equating the result with (6c), gives 

A " k = W ( a _> CM * " + M 9 It C n m ( 0 k>\) ' S n m (0 k ,9 k )] (& 


1 r k 


r t n-1 


nT = W 'a ' n s n m < 9 k ♦ Mk-^-c„"'(e k ,^) J , 


where (M k , Mq, m|) are the components of M k . 
Detailed models of the phase information 

g nx - , \ . h n m = z B n £ 


are not needed to derive theoretical forms for the expected crustal magnetic spectrum 
a 2n+4 n 

IVMI = (»+»(-) ( I <8„?) 2 ♦ <h n ™) 2 | , 9bl 

m=0 v ' 

- HW&W 2 * 44 (?c) 

implied by simple sialistical hypoiheses about the K domain dipoles. Indeed, by (9a-c), the reference 
spectrum for the crustal field produced by these sources is expected to be 

{R nx } - (n+I){ m y ^An?) 2 M X^) 2 , (1Qa) 

= ^ + J) { ml l li if, ^ An?) + (B " k B n * )] 1 ■ dOb) 

The functional dependence of the coefficients in (lOa.b) upon the domain dipole moments is given by 
(8a,b). Detailed models (9a) are, however, required to test statistical hypotheses about crustal 
magnetization by comparing expected spectra (10a,b) with observational spectra (2b) (see section 6). 

2.1 Spectrum of a Single Offset Dipole 

The magnetic spectum of a single point dipole M^r^) is 


6 



( 11 ) 


a 2n+4 n 2 2 

R nx« - (" +D(-) Z (A n £) + (Bj) . 

r m=0 


By (8a,b), the sum over order m in (1 1) is 


n 

Z 

m=0 


( A nk) 2 


+ (B 


n\\2 _ 
nk' 


n n j 2 

G n Z (Mk) 2 n 2 [P n m ]2 + G n Z (M§) 2 [ P n m ] , 

m=0 m=0 d0 k 


+ G 


n 

Z 


M. 


n - (M£ 
m=0 


m 


sin 2 ©^ 


[P n m ] 2 + G„ 

m=0 d0 k 


P m l 
n J 


( 12 ) 


where G n (r k ) _ (4jia 3 )* 2 (r k /a) 2n ' 2 and it is understood that both P n m and its derivative are evaluated at 
cos0 k . 

To cast (11) into a more illuminating function of M^, five harmonic identities (13a-e) are used to 
reduce the four sums over m on the right of (12). For arbitrary (0,<|>), the sum rule for the Schmidt 
normalized harmonics 


Z [C n m (0,<t>)] 2 + [S n m (0,<|»] 2 = X (P n m ) 2 = 1 
m=0 m=0 


(13a) 


(see, e.g., Jackson [1975, equation (3.62) with y = 0] orLangel [1987]) is repeatedly differentiated with 
respect to 0 to obtain 


I - 2[ ® p ""’ ltP "' n| - 0 


m=0 


(13b) 


n d 2 d 

X 2 [ P n m ][P n m ] + 2 [ ~ P n m ] 2 = 0. 
m=0 d0 2 d0 


(13c) 


The first sum (of radial terms) on the right of (12) is (M*) 2 n 2 by (13a); the fourth sum (of cross terms) is 
zero by (13b). 

To reduce the second sum (of colatitudinal terms) on the right of (12), recall that for either C m or 
S n m the surface Laplacian is n 

o 1 a as n m 1 a^ro 

V s 2s n m - -n(n+l)S n m = — ± s in0 — + 

sin0ae 30 sin 2 0 34,2 ’ 


2d m 


d z P 


d© 2 


n m z m 

— = [-n(n + 1) + ]P n m 


cos0 ^n™ 


sin z 0 


sin0 


d0 


so 


n m 2 


2 cos0 “Pn" 1 


Z f^ P n m ]Pn m = I !-==- - ■C-.+ lHIPn-l - — [ 


m=0 d0 2 


m=0 sin 2 0 


sin0 


d0 


-][P n m ] • 
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The last terms on the right of this expression sum to zero by (13b); substituting the remainder into (13c) 
gives 

n d 2 n 2 

* [~ P n 1 = S [n(n + 1) - — — ] [P n m ] 2 . (13d) 

m =0 d0 m=0 sin 2 0 

The sum over m of [m(P n m )/sin0] 2 is needed to reduce third sum on the right of (12) as well as to 
complete reduction of the colatitudinal terms. To obtain this sum, recall that, with y being the angle 
between the position vector to the observation point r and the position vector to the source r k 

cosy = cosOcosOjj + sin0sin0 k cos(<f> - 1 )^) , 

the full addition theorem for the Schmidt normal harmonics is 

n 

ro fo + VKfc»s n , "<M k ) - P„(cosy), 

where P n is the Legendre polynomial of degree n. Repeated operation with (sin0 k ) _1 3 gives 

" ' m2 1 3 2 P (cosy) 

^0 + - -Z 

s,n e k sin\ a^ 2 

which, when evaluated at y = 0, is the desired sum. Because 

^Pnfcosy) d 2 P n (cosy) acosy2 dP„(cosy) ^ 

" — * [ ] + L 

d^ 2 dcosy 2 ^ dcosy % 

and noting that 


1 3cosy . 

— r— = sin0(sin0 k )* 1 sin(<|> - (k.) 

sin 2 e k <*Pk 


is zero at y=0 while 


1 3 2 cosy 
sin0 2 k 34^.2 


-sinOfsinOjjJ^cosfiJ) - 4^) 


is -1 at y = 0, the desired sum reduces to 


" -m 2 dP n (cosy) 

I — - [P n m (cos0 k )] 2 

m=0 sin^G^ dcosy 


= -dP n (x)/dxl x=0 = -n(n + l)/2 . 
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In (13e), x denotes cosy and the final step follows from the Legendre equation at x = 0 and the usual 
normalization of Legendre polynomials (P n (0) = 1). 

By (13e), (13d) is equal to n(n+l)/2. The second sum on the right of (12) is thus n(n+l)(M^) 2 /2 . 
Moreover, by (12e), the third sum (of longitudinal terms) on the right of (12) is n(n+l)(Mjp 2 /2. 
Substitution of (12), as simplifed via identites (13a-e), into (11) yields 

1 r^ 2n-2 j 

R nx k W = TV ( ~ } tn 2 (n+l)[M j] 2 + - n(n + l) 2 [(M§) 2 + (M J) 2 ] ] . (14) 

(47tr*r r 2 ^ 

Equation (14) gives the spatial magnetic power spectrum on the sphere of radius r due to a point dipole of 
moment M k located inside the sphere at position r k . Note that there is no radius at which (14) becomes 
independent of degree n. Compared with the plain exponential spectrum that levels off at radius r k (R n * = 
A (r k /r)^ n+ ^), the spectrum of an offset dipole increases as a cubic polynomial of n. 

Spectrum (14) is insensitive to the north-south vs. east-west orientation of a dipole because R n (r) is 
invariant under rotations of the coordinate system. Yet a well-offset radial dipole can produce almost twice 
the mean square field (E n R n (r)) of a horizontal dipole because R n (r) does vary under coordinate 
translations. This seems intuitively obvious in light of the inverse cube law and the anisotropy whereby, at 
a fixed distance from a dipole, field intensity on its axis is twice that on its equator. The area of strong field 
intensity directly above a well-offset radial dipole clearly dominates the mean square field on the sphere. 
For a horizontal dipole of the same absolute moment at the same modest depth, there are two areas of 
strong field on the sphere; however, these areas are further from the dipole itself, have lesser intensity, and 
thus give a smaller mean square field than does the radial dipole. Of course, a geocentric magnetic dipole 
of moment m^ is purely radial and, with r k = 0, equation (14) correctly gives the single non-zero R n , Ri = 
2(p 0 m 0 /47n 3 ) 2 . ° 


2*2 Randomly Oriented Dipoles on a Spherical Shell (Spectrum Iss) 

At a particular instant in geologic time, the K domain dipoles are regarded as a sample population 
with mean square moment { M 2 ) . The coefficients appearing in (10) are given in terms of domain dipole 
moments M k by (8a,b). The positions of the dipoles are taken to be random samples from a laterally 
uniform spatial distribution. The positions are thus expected to be uncorrelated in that {r k »rj} is equal to 
( (ifc) 2 } , where the Krtfnecker delta 5^ is 1 if k = i and is 0 if k * i. 

If the dipole moments are randomly oriented, then any particular orientation is as likely as its 
opposite 


|Mk) = (Mg) = (M jp = 0, 

there is no reason to expect cross-correlated 

= {M&M*} 

there is no reason to expect cross-correlated 

(M j M j) = {M q M 

= {M&Mj) 


(15a) 

components for an individual moment 


= 0, 

(15b) 

moments 


= 0 for k * i 

(15c) 

= 0 for k * i , 

(15d) 


but the auto-correlations remain perfect 
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< M r M r> = {M§M&} = (m|m|) = {M 2 )/3. 


Equations (15a-e), summarized by {M*} = 0 and [M j k M,M = { m 2/ 3 J provide a mathematical 
statement of the random dipoles hypothesis. 

^Although M k depends upon position (8^), position is independent of M k ; therefore, in evaluating 
I A n k } and (B n k ) via (8a,b), the expectation operator passes right through the harmonic functions of 
(0 k 4>k). For random dipoles the expectation values of the coefficients are thus zero: (A ) and { B Pi 
are zero by (15a) and |(8a-b)|, so {g n ^| and (h n ^| are zero by {(9a)}. It follows that the expected field 

components are also zero. Moreover, (15c-d) and (8a-b) imply that contributing coefficients are 
uncorrelated 

{A nkA n ?> = {B n <£B n m} = 0 fork*i. (]6) 

The expected field vector from random dipoles is zero, but the expected spatial power spectrum, 

hence the expected field intensity, is not. Substitution of (16) into (10b) yields the expected reference 
spectrum 

n K 

{R nx> (n + 1) { I [ I (A n j?) 2 + (B n £)2]j. 

m=0 k=l u/ ' 

Tlie order of summations in (17) is reversible, so the expected magnetic spectrum from K random dipoles is 

the sum of the expected spectrum from each dipole. If all K dipoles are at radius r_, then (14) and (17) 
imply * 


liss _ 


r x 2n-2 


^y ( 7 ) [nVl)|M r 2 l + i 2 nO+ D 2 IMe 2 + M^H . 


With (15c), (18) simpifies to 


3(4na 3 ) 2 


[n (n+1) (2n+l)] ( — ) 


r x 2n-2 


TTiis is the reference spectrum expected from K random dipoles of mean square moment {M 2 } on the shell 
of radius r x . With amplitude defined via aJ s (2/3)K{M 2 }/(4jca 3 ) 2 (19) is rewritten as 

^nx^ SS = A x n (n + 1/2) (n + 1) (r x / a) 2n ' 2 (20a) 

and the expected magnetic spectrum above the spherical shell is just 

< R nx«> ISS = {R„ x ) Iss (a/r) 2 ^. (2Qb) 

There is no radius at which spectrum (20) becomes independent of degree n. Compared with an 
exponential spectrum that levels off at r x , the spectrum expected from a shell of random dipoles increases 
as a cubic polynomial of spherical harmonic degree n. 

Although n is a positive definite integer, the derivative of (20b) with respect to n vanishes when 
(2n 3 + 3n 2 + n) ln(r x /r) = -(3n 2 + 3n + 1/2) ; 


1 1 



therefore, spectrum (20b) has finite extrema if and only if r x < r. To the extent that Earth’s magnetic crust 
resembles a randomly magnetized shell of radius r x < a, its reference spectrum (2c) should resemble (20a). 
If the positive definite shell depth d = a - r x is small in that d « a, then reference spectrum (20a) has a 
maximum or ’peak’ at high spherical harmonic degree N xp « -3/[21n(r x /a)] - 1.5(r x /d) With a = 6371 km, 
shell depths of 10 km, 20 km, or 40 km respectively give N xp of 954, 477, or 238. 

The geometric attenuation factor makes the degree of the peak of spectrum (20b) decrease as radius r 
increase. For example, at satellite altitude 450 km above, and a shell 10 km below, the reference sphere, 
the peak is at harmonic degree 21. 


2,3 Random Polarity Field- Aligned Dipoles on a Spherical Shell (Spectrum Uss) 

Thermoremanent magnetization (TOM) of volcanic rock cooling in the aperiodically reversing 
geocentric axial dipole field from the core suggests non-randomly oriented domain dipoles. The crustal 
magnetic spectrum appropriate to field-aligned domain dipoles derived in this section is for domain dipoles 
that are either parallel or anti-parallel to the present axial dipole field with equal liklihood. Subsequent 
reorientation of domain dipoles by crustal deformation, metamorphic remagnetization, sea-floor spreading 
and continental drift, and/or true polar wander may lead to a crustal spectrum more akin to that derived in 
the previous section; however, the two spectral forms turn out to be remarkably similar. 

If the domain dipoles are not randomly oriented, but are either parallel or anti-parallel to a geocentric 
axial dipole field, then 


M r “ 2M^cos0 k , = M*sin0 k , m| = 0. 


( 21 ) 


and M can be either positive or negative. With no preferred polarity, each M £ is a random sample from a 
population with zero mean value, (M = 0 and 


(M*} = |M{$) = |M|) =0. ( 22a) 

For example, the PDF for M £ could be symmetric due to remanent acquistion processes that are 
independent of core dipole polarity which, in tum, is independent of absolute core dipole moment. With no 
reason to expect cross-correlated polarites, { M £ ) is zero for k*i,so 

= { 2M ^ cos0 k sin0 j } 

= cos0 k sin0 i = 0 for k ^ i 


and the expected cross-correlations are zero 


= 0 fork*i 
= 0 for k * i . 

The auto-correlations remain perfect, but now depend on colatitude 
(M*^ M^l = 4M Q 2 cos 2 0 k 
= M Q 2 sin 2 0 k 
= 2M Q 2 cos0 k sin0 k 


(22b) 

(22c) 

(22d) 

(22e) 

(22f) 
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where M Q 2 is the ensemble average {(M k ) 2 }. The latter is independent of position, so the sum of (22d) 
plus (22e) can be readily integrated over the sphere to obtain <{M k *M k }> = 2M Q 2 . 

For random polarity field-aligned dipoles the expectation values of the coefficients are again zero: 
{Ank* I ^ f^nk I 216 zero (22a) and { (8a-b) } , so {g n ™} and {h^j are zero by {(9a)}. The expected 
field components are again zero. Equations (22b-c) and (8a-b) imply that individual coefficients are 
uncorrelated (see (16)), so the expected spatial magnetic power spectrum is again given by (17). 

From (8a,b), (17), (22) and the supposition that all domain dipoles are on the sphere of radius r one 
obtains the expected reference spectrum 




IIss 


G n ( f x ) {(M 0 k ) 2 } 


X X n 2 (n+l)[2cos0 k P n n1 ] 2 + (n+ l^sin^ P m ] 2 

k=l m=0 d0 k 


or, by identites (13a-e), 


+ 


n(n + 1) cos0 k sin0j { [P n 


m p mi 

dOk n 


(23a) 


(^nxl^ SS ~ G n (r x ) M 0 2 X n 2 (n+l)[4cos 2 0 k ] + n(n + 1 ) 2 [sin 2 0j c ]/2 . (23b) 

k=l 

With domain positions amounting to random samples from a uniform distribution over the spherical shell, 
the sum in (23b) over all K positions is equivalent to the average over the sphere multiplied by K 

K 

X [n(4cos 2 0 k ) + (n + l)(sin 2 B k )/2] 
k=l 

K ^ 

= — J [4ncos 2 0 k + (n+l)(sin 2 0 k )/2] sin0 k d0 k 


- (K/3)[4n + (n+1)] = (5K/3)(n + 1/5) ; 
therefore, 


{ R nx } IIss = G n (r x )(5KM 0 2 /3)n(n+l)(n+l/5) . 

With M q 2 = {M 2 }/2 and amplitude a} 1 = 5K{M 2 )/[6(47ta 3 ) 2 ], (23c) is rewritten 

{ R nx } IIss = A”n(n+ 1/5) (n + 1) (r x / a) 2n ' 2 
and the expected magnetic spectrum above the shell is 
{ R nx (r)) IIss (r) = (R nx} nss (a/r) 2n+4 


(23c) 


(24a) 


(24b) 


Like (20a), reference spectrum (24a) increases as a cubic polynomial of n. Like (20b), there is no radius at 
which spectrum (24b) becomes independent of degree n. 

The ratio of the expected spectrum for random polarity field-aligned dipoles (24) to that for randomly 
oriented dipoles (20) is 


12 


ITT' 



t R nx> I1SS/ { R nx> ISS = ( A x n / A x I )K n+1 / 5 )/(n+l/2)]. (25a) 

which increases by but 2.5% as n is increased from 12 to infinity. This should make it difficult to 
distinguish between these two source distributions by examing observational values of R n alone. 

Nonetheless, if the sum of the squared domain dipole moments are the same for these two different 
dipole distributions, then K { M 2 ) in A x ! is the same K { M 2 } as in A x n and 

l R nxk S M R roJ SS ) = (5/4)[(n + l/5)/(n + 1/2)] . (25b) 

All else being equal, the high degree multipole power from random polarity field-aligned dipoles is thus 
about 25% more than for randomly oriented dipoles. This is attributed to the nearly radial orientation and 
typically strong moments (2M 2 ) ^ of field-alligned domain dipoles located near the poles. These 
contribute more to multipole powers R n than do randomly oriented dipoles of ordinary strength {M 2 } 1 / 2 . 
The surplus power from the polar regions evidently exceeds the power deficit due to the nearly horizontal 
orientation and typically weak moments {M 2 /2} 1/2 of the dipoles located near the equator. 

3. Theoretical Magnetic Spectra for a Spherical Annulus 
Earth s magnetic crust is not a spherical shell, so the spectra derived above should be corrected for 
effects of non-zero crustal thickness and ellipticity. To see the effect of such corrections while staying 
close to the geologic null hypothesis, the position of each dipole is treated as a random sample from a 
spatially uniform distribution within a magnetic crust of uniform thickness 2d. 


3.1 Random Dipoles in a Spherical Annulus (Spectrum Isa) 

As noted in section 2.2, the expected magnetic spectrum from K random dipoles is the sum of the 
expected spectrum from each dipole. This remains true if such dipoles are distributed in a spherical annulus 
of inner radius r x - d and outer radius r x + d. So (18) is here replaced by 


< R nx> 


Isa = 


{M 2 } 

3(4jra 3 ) 2 


K r k 2n ' 2 

E ( — ) n (n + 1) (2n + 1) . 

k=l a 


(26) 


The sum over k is equivalent to K times the volumetric mean, so 


^ R nx) 


Isa _ 


4jcK(M 2 } r x +d 

— n (n + 1) (2n +1)1 

3(4na 3 )X 

K{M 2 } r x +d 2n+1 

— — n (n + 1) [( ) 

3(4na 3 ) 2 X a 


r ’ 2n-2 

( — ) r’ 2 dr’ 

a 


r x -d 2n+1 


a 


(27a) 


(27b) 


where the volume of the annulus is X = (4n/3)((r x +d) 3 - (r x -d) 3 ]. As anticipated, (27b) reduces to (20a) in 
the limit as (d/r x ) approaches zero. The alternate form 


< R nx> 


Isa _ 


K{M 2 ) r x +d 2n+1 

— n (n + 1) ( ) [l 

3(4na 3 )X a 


r x -d 2n+l 

r> 

r x +d 


(27c) 


13 



shows that if the top of the layer is at the reference radius, then expected reference spectrum (27) at very 
high degree is proportional to n 2 rather than n 3 . 


With r x + d = a, A* a = A 1 (4;ia 3 /2X), and 


< R nx! ISa = Aj a n(n+l)[l - ( 


r x -d 2n+1 


-) ] 


t a-9r! 2n+l 

= A* a n(n + 1)[1 - ( ) ], 

a 

the expected magnetic spectrum above such an annulus is 
« R nx( r )) ISa - (R nx ) Isa (a/r) 2n+4 . 


(28a) 


(28b) 


Reference spectrum (28a) is equivalent to that obtained by Jackson [1990, equation (21)] using a more 

elegant method. Note that for a solid ball of random dipoles, the expected magnetic spectrum is directly 
proportional to n(n + l)(a/r) 2n+4 . 

The ratio of thick crustal spectrum (28a) to thin shell spectrum (20a) is the thickness factor 


2n+l r -d 2n+l 

T n = {Rnxl ^dRnxl 155 = (4Tcr x 3 /X)(2n-H 1 )' 1 [(— ) - (— ) ]. ' (29) 

r x r x 

For small relative thickness 6 = 2d/r x « 1 and lower degrees n «8 _1 , T„ is close to unity and (27) can be 
approximated by (20). Such approximation fails at higher degrees near or exceeding 8' 1 . For example, the 
first, second and third columns of Table 1 list n, 360°/n, and the values of T n for a 40 km thick crust 
obtained with a = 6371 km, r x = 6351 km, and 8 * 1/159. The thickness correction amounts to 10% at 
degree 120. Therefore, a 40 km crustal thickness might be distinguished from an equivalent source shell 
depth of 20 km - provided the R n could be determined with 10% accuracy for degrees through 120. Fitting 
expected spectrum (28a) to similarly uncertain R n values from a degree 60 model of satellite data is, 
however, unlikely to give a very much better estimate of magnetic crustal thickness 2d than is fitting (20a) 
and concluding 2d = 2(a - r x ). 


3.2 Random Polarity Field-Aligned Dipoles in a Spherical Annulus (Spectrum lisa) 

For field aligned dipoles of random polarity distributed in a spherical annulus of inner radius r - d 
and outer radius r x +d, the analog of (27) is x 


(R 


nx 


jllsa _ 


5(4ttK)M 0 2 

3(4rca 3 ) 2 X 


r x +d r > 2n-2 

n (n + 1/5) (n + 1) / (— ) r >2 ^ 


5K{M 2 } 

6(4ra 3 )X 

5K{M 2 } 

6(4tta 3 )X 


n (n + l/5)(n + 1) 
(2n + 1) 


n (n + 1) 


(n + 1/5) 
(2n + 1) 


r x -i-d 2n+ l 

K ) 

a 

r x +d 2n " 2 

( ) 

a 


r x -d 2n+1 

- ( ) 1 

a 


[1 


r x -d 2n+1 

(— > ]. 

r x +d 


(30a) 


(30b) 


(30c) 


Table 1: Sample Thickness and EUipticity Correction Factors for 


n 

type I and type II dipoles (T n 
360% T n (40km) Q n 

based on a = 637 1 .000) 

E 1 E n 

n c n 

TE I 
in n 

TE H 

12 

30 

1.00091 

0.9737 

1.00023 

0.9857 

1.00116 

0.9860 

18 

20 

1.00208 

0.9611 

1.00056 

0.9786 

1.00267 

0.9801 

24 

15 

1.00373 

0.9487 

1.00103 

0.9865 

1.00480 

0.9748 

30 

12 

1.00586 

0.9366 

1.00164 

0.9718 

1.00755 

0.9702 

36 

10 

1.00847 

0.9248 

1.00239 

0.9586 

1.0109 

0.9662 

48 

7.5 

1.0151 

0.9021 

1.00430 

0.9462 

1.0196 

0.9600 

60 

6 

1.0238 

0.8803 

1.00676 

0.9344 

1.0307 

0.9561 

72 

5 

1.0344 

0.8595 

1.00976 

0.9234 

1.0445 

0.9546 

90 

4 

1.0541 

0.8301 

1.0153 

0.9081 

1.0702 

0.9567 

120 

3 

1.0975 

0.7853 

1.0271 

0.8857 

1.1127 

0.9715 

180 

2 

1.228 

0.7090 

1.0605 

0.8518 

1.302 

1.0449 

240 

1.5 

1.425 

0.6470 

1.107 

0.8310 

1.577 

1.183 

300 

1.2 

1.708 

0.5960 

1.166 

0.8217 

1.991 

1.402 

360 

1 

2.102 

0.5536 

1.239 

0.8230 

2.602 

1.728 

480 

0.75 

3.383 

0.4876 

1.427 

0.8549 

4.819 

2.886 

600 

0.6 

5.767 

0.4390 

1.680 

0.9240 

9.667 

5.315 

720 

0.5 

10.23 

0.4018 

2.011 

1.0313 

20.51 

10.51 

900 

0.4 

25.38 

0.3598 

2.694 

1.2728 

68.08 

32.16 
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where X is again the volume of the annulus. The form (30c) again shows that, when the top of the layer is 
at the reference radius, the expected reference spectrum at very high degree is proportional to n 2 With r + 
d = a. A* 13 s A x u (4;ta 3 /2X), x 


(R, 


nx’ 


lisa _ 


A x IIa n(n+ 1) 


(n + 1/5) 
(n + 1/2) 


r -d 2n+ * 

[1 - ( ) ] , 


(31a) 


and the expected magnetic spectrum above such an annulus is 

{R n x<r» IIsa = {R nx ) IIsa (a^) 2 n^. (3 , b) 

The difference between thick crust spectrum (31) and thin shell spectrum (24) again remains small until 
harmonic degree n approaches the reciprocal relative thickness (n * r x /2d). Indeed, the ratio of these 
spectra for type II dipoles gives the same crustal thickness factor T n as for type I dipoles (see (29) and third 
column of Table 1). This is becasue any systematic change in the absolute moment of type II dipoles with 
depth, perhaps due to the increase in the strength of the alligning field, has been omitted for simplicity. 


4. Theoretical Magnetic Spectra for a Thin Ellipsoidal Shell 

Consider K domain dipoles on an ellipsoid of revolution about the polar axis r e (0). The radius of the 
k n dipole depends upon colatitude Qj. 

MM 2 = ^ x 2 (1 - e 2 sin 2 0 k )' 1 = a x 2 (l + Y 2 cos 2 e k r 1 , (32a) 

where a x is the semi-major axis, b x is the semi-minor axis, e is the (numerical) eccentricity e 2 * 1 - 
(b x /a x) ’ and y is here called the "coeccentricity" y 2 = (a x /b x ) 2 - 1 = e 2 /(l-e 2 ). The eccentricity of 
oblate spheroidal shell (32a) is here taken to be that of the reference ellipsoid. For the latter, Nerem et al 

[1994] give a e = 6378.137 km and b e = 6356.7523 km; therefore, e 2 = 6.694384X10- 3 . 

The volume of the oblate spheroid is 

n 

V e = (2nb x 3 /3) J q (1 - e 2 sin 2 0)-3/2 s ined0 = 4rta x 2 b x /3 . ( 32 b) 

For a e and b e above, (32b) gives a volume equal to that of a sphere of radius 6371.001 km. The differential 
surface area element of the ellipsoidal shell fc 


dA e = 


e 4 sin 2 0cos 2 0 ^ ,, 

+ “ — T~r] [r e (0)] 2 sin0d0d<|» 

(1 - e 2 sin 2 0) 2 


= 0 - e 2 sin 2 er 1 [l - e 2 (2 - e 2 )sin 2 0] 1 ( 2 [r e (0)] 2 sin0d0d<|> 

= (1 + y^os 2 ©)-^ + a 2 cos 2 0] 1 ( 2 [r e (0)] 2 sin0d0d<t> 

= ( 1 + y^x 2 )- 1 [ 1 + a\ 2 \ [r e (x)] 2 dxd(j> 

where a ■ (y 2 + e )/(l - e 2 ) = 2y 2 + and x a cos0. The area of the oblate spheroid is 


^e “ 


2k J [1 + 
0 


e 4 sin 2 0 cos 2 0 ^ 2 
( 1 - e 2 sin 2 0) 2 


[r e (0)] 2 sin0d0 


(32c) 



= 2rca x 2 


+ 



1 +e 


1 -e 


■) 


(32d) 


For a e and b e cited above, the radius of the sphere with area A e is 637 1.007 km. 

4.1 Random Dipoles on an Ellipsoidal Shell (Spectrum Ies) 

Statistical properties of random dipoles are still described by (15). The reference sphere of radius a 
does not enclose those sources at b x < r k < a; yet R„ x k (a x ) is merely R nx k (a)(a/a x ) 2n+4 by (14). 
Combining (14), (17), and (15) now gives 

i {M 2 } K 2n ' 2 

< R nx) = — tn (n + l)(2n + 1)] I (— ) , {33) 

3(4jta 3 ) 2 k=l a 

where r k depends upon Oj. by (32a). 

The sum over k is equivalent to K times the average over the ellipsoidal shell, so the reference 
spectrum is 


j» lies 2K{M 2 )47ta 2 
l R nxl ~ TT n (n + 1/2) (n + 1) I n ; 

3(4rca 3 ) 2 A p 


(34a) 


the expected mean square field per harmonic degree n averaged over a sphere above the shell is 

UW r » IeS = (R nx } Ies (a/r) 2n+4 ; 

and, in terms of x = cosG, the integral in (34a) is 


(34b) 


1 r e (x) 2n (1 + a 2 x 2 )^ 


I„ * I l— 1 [• 

0 3 (l+^x 2 ) 


-] dx . 


(34c) 


With [r(x)] 2 = a x 2 /(l + y^x 2 ) from (32a), 

\ 2n 1 , t n 

i„ - <-) i d*^.v VaV®* 
a 0 

* (a x /a) 2n Q n . 


(35) 


Integrals Q n were evaluated numerically using Simpson’s Rule. An analytic recursion relation that 
accurately approximates the Q n is derived in Appendix B. 

To correct the spherical shell reference spectrum (19) for ellipticity, simply multiply it by the 
ellipticity factor 


E n * f R nxl IeS /( R nx> ISS = (aA x ) 2n ‘ 2 (4 ra 2 /A £ ) Q n (a x /a) 2n = (a^) 2 " (4rrr x 2 /A e ) Q n . (36) 

If r x * s chosen such that 4irr x 2 = A e , then E n ^ is simply Q n (a x /r x ) 2n and the increase in (a x /r x ) 2n with n 
partly compensates the decrease in Q n . 
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Table I shows values of E n J computed using the semi-major and semi-minor axes (a p , b ) cited 
above and the radius r x of the sphere with area equal to that of the corresponding reference ellipsoid. The 

columns of Table 1 give: n, 360% T O f I - n /, /, \2n ru:\ A .. ,. 

, ' * re %• % “ "n^W front (36) and other terms discussed in 

o owing sections. E n exceeds unity because elevating random dipoles on the equatorial bulge adds more 

to |R nx } above the shell, hence to E n , than depressing dipoles near the poles subtracts from it. As 

anticpated m light of section 3.2, the corrections are appreciable (+17% or more) at harmonic degrees near 

or exceeding the reciprocal flattening (n > 1/f = a x /(a x - b x ) * 298). 


4.2 Random Polarity Field-Alligned Dipoles on an Ellipsoidal Shell (Spectrum lies) 

Field alligned domain dipole moments M k should vary systematically over the oblate spheroidal shell 
(32a), not only because an alligning geocentric axial dipole field varies with colatitude, but because the 
(a/r) attentuation further weakens such a field on the equatorial bulge and strengthens it on the flattened 
po es. The latter indicates modification of (21); however, omitted variations in the alligning field (dipole 
tilt, non-dipole core fields, and temporal changes) are already larger than this effect. Equations (21) and 

(22a-f) are therefore retained and, with (8a,b) and (10b), the expected spectrum remains the sum of K 
contributing spectra (17). 

As in section 4.1, ellipticity prevents attenuation factor G„(r k ) = (4*a 3 r 2 [r (0, )/ a ] 2 "- 2 from beins 
replaced with G n (r x ). So (23b) is replaced with 


(R )Hes = 

\ JY nx ) 


Mp 2 

(4ita 3 ) 2 


K r k 
n (n+1) X (— ) 
k=l a 


2n-2 


[4ncos 2 6 k + (n + l)(sin2e k )/2)] . 


(37a) 


With (32a) for r k = r e (0 k ) = r e (x k ) on the ellipsoid. 


Mq 2 


\ 2n ' 2 K 


|Rn,l|neS * ^2 n< " +1) <r f ) k !, [4, “ t2 * ( n + 'X 1 * X]t 2 )/2K1 + At 2 ) 1 *" . 


(37b) 


Again, the sum is replaced with K times the average over the shell, which amounts 

differential area element (32c) and dividing by the area of the ellipso idal shell. The 
spectrum 


to integrating over 
resulting reference 


{r (lies _ 


'nx> 


2tta 2 K{M 2 ) 

(4na 3 ) 2 A p 


n(n+l)J n , 


with 

a x 2n 1 

J n * ( — ) J [4nx 2 + (n+lKl-xVKl + ^xVn-id-aVj^dx 
a q 


5n + 1 a x 2n 

— ( T } Yn 


(38a) 


(38b) 


becomes physically meaningful above the shell (r > a x ) where 
{ R nx ( r )) IIes = (R nx > IIes (a/r) 2n+4 . 
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Comparison of (38a-c) with (23c) shows the ellipticity correction factor for random polarity field 
alligned dipoles to be 

E n n = ( R nx) IIeS /f R nxl IISS = (V r X ) 2 " ( 4rtr x 2 /V Y n • (39) 

Table 3 shows values of E n ^ computed using numerically integrated Y^, the semi-major and semi-minor 
axes given above, and r x = (A^n)^ 2 . Unlike E n ^, correction factor E n ^ is less than one for low through 
intermediate degrees (n < 660). This is because ellipticity places the stronger, nearly radially oriented 
domain dipoles near the flattened poles (see section 2.3) and thus further from a source-containing sphere 
on which the spectrum is evaluated. The corrections are appreciable at harmonic degrees near or exceeding 
the reciprocal flattening 298 (e.g., -18% at n =300). Correction factor E n ^ begins to increase with n above 
this degree and exceeds unity at degree 720. This reflects the dominance of the geometric factor (a x /r x ) 2n 
at sufficiently high degree: elevation of weaker horizontal domain dipoles on the equatorial bulge 
ultimately adds more to very high degree multipole power, hence to E n n , than depressing stronger radial 
dipoles on the flattened polar caps can subtract from it. 

5. Theoretical Magnetic Spectra for an Ellipsoidal Annulus 
The expected spatial magnetic power spectrum for an ellipsoidal annulus of either randomly oriented 
(type I) or random polarity field-alligned dipoles (type II) is here obtained for an oblate layer of invariant 
eccentricity e and thus slightly variable thickness. In particular, at the top of the shell 

[ r e(®k^ 2 = b x 2 (l-e 2 sin 2 0 k r 1 = a x 2 (1 + y^cos^)' 1 (40a) 

and at the base of the shell of equatorial thickness 2d = a x - 

[ r e( 0 k^ 2 = b* 2 (1 - e 2 sin 2 0 k )' 1 = (a x -2d) 2 (l + y 2 cos 2 e k )' 1 , (40b) 

where b* = (^ - 2d)(b x /a x ). The volume of this annulus is 

x ea = (47t/G)[a x 2 b x - (ax - 2d) 3 (b x /a x )] = 47ta x b x (2d)[l + 2d/a x - (2d 2 /a x ) 2 /3] . (40c) 


5*1 Random Dipoles in an Ellipsoidal Annulus (Spectrum lea) 

The statistical properties of random dipoles are still described via (15). Combining (14), (17), and 
(15) again gives (33), but the range of r k is given by (40a-b), For a spatially uniform distribution of domain 
dipoles, the sum over K dipoles is replaced by K times the volumetric mean. The differential volume of 
integration is drdA e , with dA e given by (32c); integration over $ yields 


( R nx> 


lea _ 


2jta 2 K(M 2 ) 1 

— n (n + 1) (2n + 1) J 

3(4na 3 ) 2 X ea -1 



r 2n (l + a 2 x 2 )^ 

(_) — drdx. 

3 (l+^x 2 ) 


(41) 


Evaluation of the radial integral at the 0 dependent limits (40a, b) yields 


t KIM 2 ! a x 2n+l a„-2d 2n+l 

{ R nx* Iea = “ n(n +!)(—) [1 - ( ) ]W (42a) 

3(4ra 3 )X ea 3 a x 
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where 


' 1+oV 


W n e J [1 +7^x^]' n ‘* [ 

0 1 + y^ 2 


J^dx 


(42b) 


The ratio of spectrum (42a) for random dipoles in ellipsoidal annulus [a x -2d, a x ] to spectrum (27b) for 
random dipoles in spherical annulus [r x -d, r x +d] gives an ellipsoidality factor 


< R nx> Iea (a x ) 2n+1 ■ (a -2d) 2n+1 

= WX ea ) 


< R nx> 


Isa 


(r +d) 2n+1 - (r x -d) 2n+1 


w„ 


(43) 


Comparison of the expected reference spectrum for random dipoles in an ellipsoidal annulus (42) 
with that for random dipoles on a spherical shell (20) leads to the thick ellipsoid factor for type I dipoles, 

TEn 1 * (4m- x 3 /X ea )( 2n + l)- 1 (a x /r x ) 2n+1 ( 1 - [(a x -2d)/a x ] 2n+1 ) W„. (44) 

The seventh column of Table 1 lists TE n * computed using a x = 6378.137, b x = 6356.7523 km, and d = 20 
km as before; however, the radius of the sphere with area equal to that of the ellipsoid with semi-major axis 
a x - d is used for r x = 6351.03 km. As anticipated, TE^ differs but slightly from the product T n E n *. 

S.2 Random Polarity Field-Alligned Dipoles in an Ellipsoidal Annulus (Spectrum Ilea) 

To obtain this final sample spectrum, the sum over K random polarity field alligned domain dipoles 
in (37a) is again replaced by the volumetric mean over the ellipsoidal annulus. The differential volume 
element is drdA e . With dA e from (32c), evaluation of the <j> integral gives 


t R nx) 


27ta 2 KM 2 
lies ° 


(4rea 3 ) 2 X 


1 r e r 2n 
n (n + 1) j J (— ) x 

-1 r" a 

ea 1 r e 


[4nx 2 + (n + I)(l - x 2 )/2] £ +oc2x2 > 1/2 ^ ^ 

(1 +y 2 x 2 ) 

Evaluation of the radial integral at the 0 dependent limits (40a, b) yields 


(45) 


< R nx> 

where 
LL s 


Ilea _ 


KM„J 


(4jta 3 ) 2n + 1 a 


n(n+l) ^ 2n+ l a x' 2< ^ 2n+l 

(-) [!-(—) ]U n 


3 1 4nx 2 + (n+l)(l-x 2 )/2 (1 + a 2 x 2)l/2 

• [ ] dx . 


5n+1 0 (1 +y 2 x 2 ) n+1 (l + fx 2 ) 


(46a) 


(46b) 


Comparison of the expected reference spectrum for random polarity, field-alligned dipoles in an 
ellipsoidal annulus with that for such dipoles on a spherical shell leads to the thick ellipsoid factor for type 
II dipoles, 


TE n 11 = (4nr x 3 /X ea ) (2n + l)' 1 (a x /r x ) 2n+1 (1 - [(a x -2d)/a x ] 2n+1 ) U„ . 


(47) 



Table 2: Elementary Spectral Forms {R nx } at Reference Radius a. 


Type 

Description 

Form 

Amplitude A in 
K(M 2 }/(4ra 3 ) 2 

Flat Spectrum 
at radius r x 

A(r x /a) 2n ’ 2 

Undefined 

Random Dipoles 



Iss 

Spherical Shell (20) 

A n(n + l/2)(n + l)(r x /a) 2n ' 2 

2/3 

Isa 

Spherical Annulus (28) 

A n(n + 1)(1 - [(r x -d)/a] 2n+1 ) 

4na 3 /3X 

Isb 

Spherical Ball (radius a) 

A n(n + 1) 

I 

Ies 

Ellipsoidal Shell (34) 

A n(n + l/2)(n + l)(a x /a) 2n Q n 

(2/3)(4na 2 /A e ) 


or (36) 

A n(n + l/2)(n + l)(r x /a) 2n - 2 E,, 1 

2/3 

lea 

Ellipsoidal Annulus (42) 

A n(n + l)[(a x /a) 2n+1 - ((a x -2d)/a) 2n+1 ] W„ 

4jta 3 /3X ea 


or 

A n(n + l/2)(n + l)(r x /a) 2n - 2 TE n I 

2/3 

Ieb 

Ellipsoidal Ball 

A n(n + l)(a x /a) 2n+1 W n 

4ra 3 /3X ea 

Random Polarity Field-AUigned Dipoles 


IIss 

Spherical Shell (24) 

A n(n + l/5)(n + l)(r x /a) 2n ' 2 

5/6 

lisa 

Spherical Annulus (30) 

A n (n + !) j; * (1 - [(r x -d)/a] 2n+1 ) 

(5/4)(4jta 3 /3X) 

Ilsb 

Spherical Ball (radius a) 

. , n + 1/5 

A n (n + 1) -i- 

n + 1/2 

5/4 

lies 

Ellipsoidal Shell (38) 

A n(n + l/5)(n + lXa^a) 2 " Y„ 

(5/6X4rea 2 /A e ) 


or 

A n(n + l/5)(n + l)(r x /a) 2n - 2 E n JI 

2/3 

Ilea 

Ellipsoidal Annulus (46) 

n+1/5 ^ 2n+l a^d 2n+l 
A " ( " +1) n + W«a> '< a > 

(5/4X4;ta 3 /3X ea ) 


or 

A n(n + l/2)(n + l)(r x /a) 2n ‘ 2 TE n n 

2/3 

Ileb 

Ellipsoidal Ball 

An(n + 1 ti (a x /a)2n+lu n 

(5/4)(4 rt a 3 /3X ea ) 


Explanatory Note: For spheres, r x denotes the radius of a shell of K domain dipoles sources with mean 
square moment { M ) ; r x -d is the radius at the base of a source layer (top at r + d = a), and X is the volume 
of die spherical annulus. Sample thickness factors T n are listed in Table 1. For ellipsoids, is the semi- 
major axis and A e is the area of the ellipsoidal shell; X„ is the volume of the constant ellipticity annulus 
spannmg [a x -2d, a x ]. Sample ellipticity factors E n ‘, E n “are given in Table 1 for 4rcr 2 = A„, as are thick 
ellipsoidal annulus factors TE* and TE n “. x e 
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The eighth column of Table 1 lists TE n n computed using a x = 6378.137, b x = 6356.7523 km, d = 20 km, 
and r x = 6351.03 km; the latter is the radius of the sphere with area equal to that of the ellipsoid with semi- 
major axis a x - d. As anticipated, TE^ 1 differs but slightly from the product of T n with E n n . 

Table 2 sumarizes the elementary theoretical forms the spatial power spectrum of a randomly 
magnetized crust derived above. 

6. Comparison of Theoretical with Observational Spectra 

We fitted some of the foregoing theoretically expected forms for the spatial magnetic power spectrum 
of Earth s crust, (R nx (r) } , to observational spectra derived via analyses of MAGSAT data, R n (r). We 
worked with the R n (r) computed from the maximum degree N = 60 model M 102 189, and the related N = 
49 model M102389, derived by Cain, Holter and Sandee [1990], At the time of our first fits (1994), these 
models were the state of the art in high degree geomagnetic field analysis. To eliminate any plausible need 
to correct for core field contributions, only R n of degrees n > 16 were fitted. Covariances for these models 
were not made available, so formally weighted fits were not possible. It was, however, possible to fit RM 
at either reference radius a = 6371 km ("at the surface”) or at radius a + 420 km ("at satellite altitude"). 
The former assigns relatively heavier, and arguably excessive, weight to the higher degree R n than does the 
latter. Before summarizing the fits, the derivation of the field models is reviewed to aid understanding of 
the origin and limitations of the observational R n . 

6.1 High Degree Models of Cain, Holter, and Sandee [1990]. 

Both M102189 and M102389 are modifications of the model M07AV6 derived by Cain et al. 
[1989a]. M07AV6 extended an N = 29 least squares model to N = 63 via numerical integration (Neumann 
or Gauss-Legendre method). Quiet (K p < 1+), near midnight, Dst corrected, hourly observatory means 
were used to obtain an improved, degree 10 secular variation (SV) model for reducing data from the 10/79 - 
6/80 MAGSAT mission to common epoch. MAGSAT data prior to 1 1/3/79 or data with K > 2 + were 
rejected, as were 4,478 of 1,330,285 observations that deviated by more than 100 nT from N P = 29 model 
M070284. Residuals relative to that model were further corrected for external field effects using (i) a Dst 
correction for ring current effects and (ii) long wavelength corrections to separated dawn and dusk data as a 
function of dip latitude for ionospheric and remaining magnetospheric effects. 

The corrected residuals were sorted into 3° longitudinal blocks centered about the roots of the degree 
64 Legendre polynomials and block mean values of AB r computed. Empty blocks (poleward of MAGSAT 
ground tracks) were filled by linear interpolation over the pole. Gauss coefficients were computed from the 
augmented block mean AB r by the usual quadrature in <)> and Gauss-Legendre quadrature in cos6. 
Coefficients up to m = 59 and n = 63 were retained; however, coefficients above degree 50 appeared to 
represent noise [Cain et al., 1989a; 1990], As noted in section 1.5, Cain et al. [1989b] interpreted the R_(r) 
spectrum from model M07AV6 in terms of a 0.091 nT 2 noise spectrum that levels off at 420 km altitude, a 
crustal spectrum that levels off 14 km below the reference sphere a = 6371 km, and a core spectrum that 
levels off about 80 km below the core surface. 

Cain et al. [1990] reexamined the 8 second MAGSAT data set underlying M07AV6 and rejected 
22,478, typically high latitude, observations that deviate by more than 50 nT from this model. From the 
remaining vector component and scalar data, 49,998 were selected to equalize the number of data per 
10°xl0° block; 276 values were synthesized from the first 50 degrees of M07AV6 to fill in blocks 
poleward of the MAGSAT ground tracks. Some abnormalities in the distribution of the resulting 50,274 

residuals relative to M07AV6 helped motivate a higher degree least squares spherical harmonic fit to these 
field values. 



Model M102189 was obtained from these values and M07AV6 by keeping coefficients with degrees 
51-60 fixed and iterating weighted least squares adjustments (with some modifications in the weights) to 
coefficients of degrees <50 (and first degree external terms). The N = 49 model M102389 was obtained by 
a single minor adjustment to a similarly truncated version of M102189. Values of R n for both models were 
generously provided by J. Cain [personal communication, 1994]. 

The coefficients of degrees greater than 50, being those of M07AV6, result from integrating block 
mean values. Two possible disadvantages to the use of block mean values are noted in passing. First, the 
calculation of a block mean value filters out information about horizontally narrow scale field structure 
within the block. Such horizontal averaging should tend to make high degree R n , as computed from a 
model fitted to block mean satellite data, underestimate true values. Second, the calculation of a block 
mean value from data acquired in even a thin annulus sampled by a satellite filters out information about 
the variation of the field with radius in each block. Such information could help establish the degree of 
contributing fields that fall off as (a/r) n+2 for large n. For example, R 50 (6371 + 400 km) is 2.15 times 
**50^71 + 4^0 km). Unlike the horizontal averaging, such radial averaging should tend to make high 
degree R n overestimate true values. For potential fields, the two effects cancel when the block is replaced 
with a small spherical ball; yet orbit decay may still cause non-uniform sampling, notably over-sampling of 
higher altitudes and a tendency for computed R n to undersdmate true values. 

6.2 Thin Shell Reference Spectra With No Noise Correction. 

The reference spectrum R n (a) computed from model M102389 [Cain et al„ 1990] is shown in Figure 
1 as the intersections of short solid line segments. Following Langel & Estes [1982], the spectrum appears 
dominated by the dipole; a broad scale non-dipole field of deep internal (core) origin for degrees 2 through 
12; and fields of crustal origin for degrees 16 and higher. Degrees 13, 14, and 15 are considered 
transitional. No special significance is attached to the behavior at degrees 2-3 or 8-9 [Voorhies & Conrad, 
1996]. The tendency of R n increase with n above 16 is qualitatively consistent with a crustal origin; 
indeed, theoretical crustal spectra (20) and (24) show {R nx ) increasing as n^. The tendency for R n to 

increase with n is also consistent with the satellite altitude noise level discussed by Cain et al [1989a b 
1990], 

Plain Exponential Fit. Piecewise exponential representations of a reference spectrum set 
ln[R„ 1 = uq + ajn ( 48 a) 

for some range of n. By analogy with (5), writing this as 
ln[R n *] = In A* + (2n+4)ln(a/r*) 

confirms that such representations level off at radius r* = exp(aj/2) and that A* = expfag - 41n(a/r*)]. 
Although no geophysical significance is here attached to r* or A*, the fit of R n * to observational R does 
provide background for the study of other spectral forms. 

Figure 1 shows the unweighted least squares fit of (48b) to the 34 values of InR^ for degrees 16-49 as 
the long solid line with positive slope. The fit appears close and amounts to 

R n * = 3.5196 (6587. 1/637 1.2) 2n+4 nT 2 . (49) 

Spectrum (49) levels off high in the ionosphere, some 216 km above the reference sphere. (This is about 
half the typical spacecraft altitude and twice the altitude of the ionospheric maximum). The sum of squared 
residuals 
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Figure 1. Spatial magnetic power spectrum on the reference sphere. R n (a) in nT 2 , as a function of 
harmonic degree n. Short solid lines segements connect observational values from model M 102389 of Cain 
et al. [1990]. Long solid line with positive slope shows plain exponential fitted to degrees 1649 (linear in 
ln(R n )). Dot-dash curve shows spectrum expected from a shell of random dipoles (20) fitted to degrees 16- 



(50) 


S 2 = ^(lnR n -lnR n *) 2 = (IndyR,/)] 2 

over the 34 harmonic degrees fitted is 1.630. Without the covariance needed to weight R n or lnR n , we are 
unable to formally test the statistical significance of the residuals. We can, however, compute the scatter 
factor f, which is the exponential of the root mean square value of ln(R n /R n *) 

f = exp[(S 2 /34)'/ 2 ] , (51) 

and is 1 .245. The R„ have been thus been fitted to within a typical factor of ( 1 .245)+ 1 . Comparison of the 
R„ from M102189 with those from M102389, again for degrees 16 through 49, gives the correpsonding 
uncertainty factor f Q = 1.0248. The residuals are therefore much too large to attribute to plausible changes 
in Field model truncation level. 

Thin Shell Fits. Figure 1 also shows the fit of expected reference spectrum (20) for randomly 
oriented dipoles on a shell to the same R n as the dot-dash curve. In particular, a linear function was least 
squares fitted to ln[R n n* 1 (n + 1/2)' 1 (n + l)' 1 ] for degrees 16-49 of model M102389. The result, 

R n ISS = 2.880x10 3 n (n + 1/2) (n + 1) (6277.1/6371.2) 2n ‘ 2 nT 2 , (52) 

gives a source shell depth of 94 km below the reference sphere. The sum of squared residuals to InR is 
now 3.206 and the R n are fitted to within a typical factor f of (1.359)+'. The n 3 dependence of the thin 
spherical shell spectrum R n ^ ss , visible as curvature in Figure 1, produces a fit which is not quite so close as 
that of plain exponential R n *. 

The similar fit of expected spectrum (24) for random polarity field alligned dipoles on a shell to the 
lnR„ of M102389, 

R n^ SS = 2.940x10 3 n (n + 1/5) (n + 1) (6276. 1/637 1.2) 2n * 2 nT 2 . ( 53 ) 

gives a source shell depth of 95 km. The S 2 of 3.221 shows that the R„ have been fitted to within a typical 
factor of (1.360) - . As noted in section 2.2, it is indeed difficult to distinguish between randomly oriented 
dipoles and random polarity field alligned dipoles from observationally accessible R n alone. 

Table 3 summarizes some spectra fitted to degrees 16-49 of M102389, and to degrees 16-60 of 
M102189, in columns listing spectral form, model, degree range, key depth, amplitude, and scatter factor. 
When either random dipoles or random polarity field alligned dipoles spectra are fitted to the higher degree 
model, the effective shell depth decreases from about 94 km to about 64 km. The shallowing of the source 
layer compensates the increase in n 3 at the higher degrees (see equations (20), (24)). For the convergent 
spectral forms, the final column of Table 3 lists the mean square crustal field on the reference sphere 
obtained by extrapolating the fit. 


100,000 

B x ims = [ r R n 'sS]l/2 
n=l 


(54) 


for spectrum Iss and the corresponding sum for spectrum IIss. 

„ A few rem arks on Figure 1 and Table 3 are offered. Firstly, the negative depth at which spectrum 
R„ levels off might be physically insignificant because the plain exponential form is unphysical and/or 
because of noise in the R n fitted. If the problem is noise, then the slightly closer fit offered by R n * is not 
physically significant. Secondly, the source shell depths from spectra R n Iss and R n " ss seem deeper than a 
mean Curie isotherm, perhaps because these physically motivated forms are inadequate, because of noise in 
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Table 3: Elementary Spectral Forms Fitted to M102189 & M102389 
with no noise adjustment 


Form 

Model 

Degrees 

Depth 

(km) 

Amplitude 

(nT 5 ) 

f 

r> ims 

X 

(nT) 

Exponential 

M102389 

16-49 

-216 

3.520 

1.245 

NA 

Exponential 

M102189 

16-60 

-208 

3.768 

1.251 

NA 

Iss 

M102389 

16-49 

94.1 

2.880xl0' 3 

1.359 

137.23 

Iss 

M102189 

16-60 

63.5 

2.188xl0' 3 

1.380 

137.56 

IIss 

M102389 

16-49 

95.1 

2.940xl0' 3 

1.360 

122.73 

IIss 

M102189 

16-60 

64.3 

2.229xl0' 3 

1.382 

123.03 



the R n fitted, or because of a loss of power from high degree R n . It turns out that a trade-off between 
amplitude A and source depth a-r x makes their coestimation very uncertain. 

The extrapolated rms crustal field listed in the final column of Table 3 is about 130 nT. This is the 
same order of magnitude as background scalar magnetic anomalies inferred from marine and aeromagnetic 
surveys. Such anomalies mainly measure the crustal field component parallel to the main field, so the total 
crustal field intensity could be larger than the scalar anomaly (by a factor of 3^ in an isotropic case); 
however, the rms crustal field ought not be less than the rms scalar anomaly. 

To check this 1 30 nT rms crustal field prediction, we analyzed the "observatory biases" derived with 
many main geomagnetic field models by R. Langel and others at GSFC. The idea is that each observatory 
bias approximates the crustal magnetic field at the observatory. R. Baldwin generously provided 10 sets of 
observatory biases obtained in the course of making 10 somewhat different geomagnetic field models. For 
each model, each observatory bias (AB = (AX~ + AY^ + AZ~) was calculated and used to compute the 
rms bias. The ten rms biases average to 488 + 106 nT. A very few observatories had very large biases, 
perhaps due to data errors, site instability, or location on exceptionally magnetic volcanic edifices. 
Rejecting the few observatory biases over 1500 nT led to 10 rms biases that average to 360 + 45 nT. The 
model using the most observatories and longest time span (1900-1995) gave the smallest rms observatory 
bias of 274 nT. All these values exceed the 130 nT prediction. The sigmas of aeromagnetic (173 nT), 
repeat (263 nT), and survey (373 nT) data relative to model GSFC(S95-sc) [Sabaka et al., 1997] also 
exceed the 130 nT prediction. It can thus be argued that the 130 nT prediction for rms crustal field is too 
low, perhaps by a factor of two or three. This might be due to overestimation of source depth and, in turn, a 
loss of power from the higher degree R n fitted. 

63 Remark on Coestimation of Amplitude, Depth and Thickness 

Scatter factors for thin shell spectra Iss and IIss to the observational spectrum (about (1.36)T^) exceed 
the thickness and ellpiticity correction factors in Table 1 (which are nearer unity than ( 1 .07)^ for degrees 
less than 60). This suggests that allowance for thickness and ellipticity ought not appreciably alter the 
apparent source depth. However, an attempt was made to estimate source layer thickness along with 
amplitude and depth; this leads to a non-linear inverse problem. 

We naively used iterative linearized least squares to fit the spectrum for a thick annulus of random 
dipoles to observational R n . The method failed to converge due to an ill-conditioned, almost singular 
information matrix. As is well known from forward modeling of magnetic anomalies, there is a trade-off 
between magnetization and thickness of a source layer. In the present context, note that spectrum (27b) can 
be written 

, To r x 2n+l 2n+l 2n+l 

l^nx I = A x n(n + l)(— ) [(l + dA x ) - (l-d/r x ) ] (55a) 

which, for sufficiently thin layers (2d « r x ) and low degrees is approximately 
t r x 2n 

f^nx I Sa “ A x a n (n + 1) (— ) (2n + l)(2d/a)] . (55b) 

The partial derivative of (55b) with respect to amplitude A x Ia is directly proportional to the partial with 
respect to d, in practice, the approximate colineanty led to a near singular information matrix. 

It was possible to fix layer thickness, estimate amplitude and depth, and then adjust the layer 
thickness to tighten the fit, but no clearly preferred solution emerged by iterating this procedure. 
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6.4 Thin Shell Spectra With Noise Correction. 

The form of the satellite altitude noise spectrum used by Cain et al. [1989b] is the plain exponential 
that levels off at 420 km. 


R„ e (r) = e [(a + 420km )/r] 2n+4 . 


Cain et al. [1989b] averaged the higher degrees (n > 50) of R„(a + 420km) to obtain the noise level e = 
0.091 nT . For a degree 60 model, the correponding rss noise would be (60e) 1/2 , or 2.3 nT. This is about 
the rss noise expected uniformly distributed attitude errors of at most 20" in the ambient field and so seems 
plausible; however, it is less than indicated by a Magsat total measurement error budget. It is not clear that 
attitude, other data errors, or modeling errors such as aliasing should cause a level spectrum (56). 

We did not expect to estimate noise level e along with the amplitude and source depth of a crustal 
spectrum due to colinearity; subsequent experiments verified this colinearity. We bypassed least squares 
coestimation in favor of mapping misfit as a function of theoretical spectrum parameters (amplitude and 
depth) at each of several noise levels. The measure of misfit computed was the standard deviation of the 
theoretical spectrum relative to the noise-adjusted unweighted observational R_(r) from degrees 16-50 of 
M102189, 


o(r) = ([32]-l X [R n (r) - {R nx ( r) ) - R^r)] 2 ) 1/2 . (57) 

n=16 

Contour plots of o at each noise level show a strong tradeoff between amplitude and depth: a shallower 
source layer can be accomodated without much change in misfit provided source amplitude is reduced. 

Figure 2 shows two examples of such o plots, both obtained from (56) at reference radius r = a = 
6371 km and e = 0.091 nT 2 . The upper panel contours the standaid deviation, as a function of amplitude A 
(nT ) and depth a-r x (km), when (R nx (r)} in (57) represents an ellipsoidal shell of random dipoles (34), 

(R nx (a)} = A n (n + l/2)(n + l)Q* n [r x /a] 2n ~ 2 . (58) 

A flaw caused Q n to exceed the Q n in Table 3 by an unimportant 0.23%. Minimum o (10.9 1 nT 2 ) is for a 
source depth of 45 km; the range of source depths giving a < 11 nT 2 is 20-55 km. This seems almost 
geophysically plausible; however, the range of source depths giving o < 12 is far larger. The tradeoff 
between amplitude and depth is clear in the upper panel of Figure 2. 

The lower panel of Figure 2 contours the standard deviation as a function of amplitude and depth, but 
with |R nx (r)} in (57) taken to be the plain exponential spectra R n * = A(r*/a) 2n+4 . The leveling depth a-r* 
is negative and a is minimum (9.57 nT 2 ) for a leveling altitude of 200 km; however, o is less than 10 nT 2 
for a broad range of altitudes between 265 km and 235 km. 

Compared with an exponential spectrum leveling off at 200 km altitude, the modest increase in a 
caused by the random dipoles spectrum (20) seems a small price to pay for a 45 km deep crustal source 

layer. Yet the misfit contoured in Figure 2 is for R„ evaluated on the reference sphere, not at satellite 
altitude. 

With no formal weights for the observational R n , misfit c(a) is dominated by the deviations at the 
higher degrees (n > 40) due to the geometric amplification upon downward continuation from satellite 
altitude. To reduce this effect, similar figures were constructed (i) for misfit measured at satellite altitude 
instead of on the reference sphere and (ii) for degrees 16-40 instead of 16-50. Many noise levels were 
sampled in both cases. The de-emphasis upon higher degree R n tends to increase the apparent depth of the 
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source layer. The main result is that a very broad range of source depths and noise levels can be 
accomodated without much change in misfit. 

For example, we summarize standard deviations relative to stallite altitude R (a + 420 km) minus 
noise level e = 0 091 nT 2 , for degrees 16-50 of M102189. The plain exponential spectrum gives o(a + 
420km) < 0. 1 nT at leveling altitudes between 120 km and 160 km; yet values < 0.2 nT 2 are found in the 
much broader range between a leveling depth of 80 km and a leveling altitude above 200 km. One local 
minimum within this range is about 20 km deep - close to the geophysically plausible 14 km leveling depth 
obtained by Cam et al. [1989b], In contrast, the ellipsoidal shell of random dipoles (57) only returns o( a + 
420km) < 0.2 nT for seemingly implausible shell depths of 60 km or more. 

To perform decisive statistical tests needed to reject various theoretical functional forms for the 
crustal magnetic spectrum, one requires reliable uncertainty estimates not only for observational R , but for 
the theoretical spectrum (e.g„ {(R n - (R n )) 2 }). In the interim, geophysically reasonable source models can 
be preferred to either elementary theoretical forms or arbitrary mathematical functions. 

6.5 Spectrum of an Equivalent Source Model of Magsat Crustal Anomalies 

Since the identification of magnetic anomalies in Magsat data, many researchers have devoted much 
effort to develop data processing techniques that help isolate the field of lithospheric origin and improve 
satellite magnetic anomaly maps for geologic interpretation (see, e.g„ Ravat et al. [1995]). One step in 
some of these techniques involves computing an equivalent source representation of reduced anomaly data. 
An equivalent source model is a network of point dipoles, typically main field-alligned and located on the 
re erence sphere, with absolute moments distributed so as to reproduce the anomaly data. In practice the 
dipoles are closely spaced, albeit fewer in number than the data, and singular value decomposition may be 
used to avoid unphysical numerical instabilites in a least squares estimate of absolute moments; 
alternatively, a damped least squares appraoch may be used. 

From one such equivalent source model with a nearly uniform dipole spacing of about 220 km, M. 
Punicker [1994, personal communication] generously provided values of the radial magnetic component 
computed at the mesh points of a 2° equiangular grid and 400 km altitude. The spherical harmonic 
coefficients and equivalent source spectrum R„<ES) through degree and order 90 were computed from these 
synthetic data by numerical integration. The spectrum peaks near degree 50, but shows ringing above 
degree 70 due to source spacing and coarseness of mesh. From apparently reliable degrees 16-55, the 
standard deviation o of simple power laws (R n (ES) ~ An9) was mapped as a function of amplitude A and 
index q. Minimum c was found for q = 1.55. Neither a level (n<>) spectrum nor a random dipoles (n 3 ) 
spectrum describes the equivalent source crustal anomaly model as well as an intermediate (n 3 / 2 ) form. 
Whether the particular q value obtained reflects correlations introduced by data reduction and equivalent 
source modeling procedures or some fundamental property of Earth’s magnetic crust is not clear. It is clear 
that a power law index less than 3 is reasonable for spatially correlated magnetizations. 

7. Summary 

The magnetic field produced by magnetization in Earth’s crust and lithosphere is distinguishable 
from the field produced by electric currents in Earth’s core becasue the spatial magnetic power spectrum of 
the crustal field differs from that of the core field. This principle has long guided anlyses of land, sea, and 
aero-magnetic survey data, which routinely subtract a broad scale field - be it a constant level, a regional 
trend, a long wavelength model, or a low degree geomagnetic reference field model - before attempting 
geologic interpretation of residual anomalies. According to the famous Langel & Estes [1982] 
interpretation of the spectrum derived from MAGSAT, the core field dominates the total spectrum, R , for 


spherical harmonic degrees n < 12 and the crustal field dominates for n > 16. Observational spectra exhibit 
two distinct regimes becasue the different radii of the two main source regions (r c « r x ) cause distinctly 
different^geometric attenuation factors for the core spectrum (r c /a) 2n+4 and the crustal spectrum 
(r x /a) n . The absence of a third spectral regime, apart from the narrow range of core-crustal overlap, 
indicates that most of Earth’s mantle is neither strongly magnetized nor carrying strong electric currents - as 
expected for hot insulating silicates. The spectral separation of core and crustal fields remains a 
geophysical distinction; it does not enable a formally unique mathematical separation of core from crustal 
field by analysis of satellite and/or near surface geomagnetic field measurements. 

The geometric attenuation factors are exponential functions of n, so a single exponential can be fairly 
closely fitted to either low, or to high, degree R n . The fit returns a radius which coarsely approximates the 
radius of either Earth’s core or Earth’s crust, respectively. This fact boosts confidence in the spectral 
separation. It does not imply that the core spectrum is flat at the top of the core, nor does it imply that the 
crustal spectrum is flat at the top of the crust. It does suggest that, after allowance for geometric 
attenuation, core and crustal spectra are sub-exponential functions of n - perhaps ratios of polynomials in n 
if not simple power laws. 

Physical models of core and of crustal sources suggest that exponential spectra should indeed be 
multiplied or divided by more slowly varying polynomial functions of harmonic degree n. For example, 
the core spectrum R nc is well described by McLeod’s rule (4e): geometrical attentuation of a spectrum that 

falls off as l/(n + 1/2) atop Earth’s core. Elementary spectra for a randomly magnetized crust are listed in 
Table 2. “ 

7.1 Two Types of Dipoles on a Shell 

Elementary theoretical forms for the crustal magnetic spectrum {R nx } derived above treat the source 
of the crustal magnetic field as an ensemble of domain dipoles. The first step derived the spatial magnetic 
power spectrum on the sphere of radius r due to a point dipole of moment M k located inside the sphere at 
position r k (14). There is no radius at which the spectrum for an offset dipole (14) becomes independent of 
degree n oHevels off. Compared with the exponential spectrum that levels off at radius r k (R n * = 
A (r k /r) ), the spectrum of an offset dipole increases as a cubic polynomial of n. 

The second step formulated the geologic null-hypothesis that the magnetic field of Earth’s crust is 
due to many domain dipoles scattered at random throughout a thin crust. Two types of dipole ensembles 
were considered. Type I, or random, dipoles are of random absolute moment and random orientation; they 
are not cross-correlated and the auto-correlation matrix is diagonal as described by (15) { M j k Mj*} = 
{M tfJSfcSj,. Type II, or random polarity field-alligned, dipoles are also of random absolute moment, but 
are oriented either parallel or anti-parallel to the modem axial dipole field with equal liklihood; they are 

also not cross-correlated, but the auto-correlation matrix is not diagonal and varies with colatitude as 
described by (22). 

The magnetic spectrum expected from many uncorrelated random dipoles is the sum of the expected 
spectra from each random dipole. The expected spectrum for type I random dipoles on a spherical shell 
(20) ts thus proportional to that of a single random dipole. There is no radius at which the expected 
magnetic spectrum from a shell of random dipoles is independent of degree n. Compared with the plain 
exponential that levels off at the shell radius r x , this spectrum increases as a cubic polynomial of n. The 
same is true of the expected spectrum for type II random polarity field alligned dipoles (24). 

The expectation spectra for type I and type II dipoles prove remarkably similar. Type II dipoles do, 
however, tend to generate more multipole power - the spectral ratio {R n II }/{R n 1 }, given by (25) for thin 
spherical shells, being about 125% for degrees greater than 12. The nearly radial orientation and typically 
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strong moments of type II dipoles in the polar caps more than makes up for the nearly horizontal orientation 
and typically weaker moments of type II dipoles in the equatorial belt. 

7.2 Crustal Thickness and Ellipticity Corrections 

When domain dipoles are scattered throughout a spherical annulus of uniform thickness 2d, instead of 
confined to a spherical shell of radius r x , the shape of the expected spectrum changes. The change amounts 
to multiplying thin shell spectra by a crustal thickness factor. The factor is quite close to unity at harmonic 
degrees less than the reciprocal relative thickness of the annulus (n < r x /2d). This is in accord with the 
common geophysical view that the length scale of a source region cannot be established until variations in 
the exterior potential field of similar scale are resolved. For a 40 km thick annulus, the correction is about 
10% at degree 120; therefore, if observational R n could be determined with 10% accuracy for degrees 
through 120, a 40 km thick crust might be distinguished from an equivalent source shell depth of 20 km. 
When fitting models to observational R n of much lesser degree and/or accuracy, a thick crust model is not 
likely to return a much better estimate of magnetic crustal thickness than simply doubling the apparent 
source depth of a thin shell model (2d = 2a - 2r x ). In practice, a near co-linearity makes coestimation of 
amplitude, depth, and thickness difficult or impossible without observational R n at sufficiently high 
degrees (n = r x /2d). 

The expected spectrum from a spherical annulus of either random dipoles (28), or random polarity 
field alligned dipoles (31), does not level off at any radius. Compared with the plain exponential that levels 
off at the top of the annulus, these spectra increase as a quadratic polynomial of n at sufficiently high 
degree (n » a/2d). This is true for a solid ball of such sources regardless of n. 

In deriving expected spectra [R nx ], the effect of oblateness is somewhat more difficult to include 
than the effect of uniform thickness. The effect amounts to multiplying the thin shell spectrum by an 
ellipticity factor. This factor is close to unity for harmonic degrees less than the reciprocal flattening (n < 
a x /(a x -b x ) = 298). For type I dipoles, the factor exceeds unity due to geometric amplification: elevating 
random dipoles on the equatorial bulge adds more to [R nx ] than depressing dipoles at high latitudes 
subtracts from it. For type II dipoles, this geometric effect also causes the ellipticity factor to exceed unity 
at sufficiently high degrees; however, there is a range of lower degrees for which the type II ellipticity 
factor is less than unity. The reason is clear: compared to a spherical shell, the nearly radially oriented 
domain dipoles located near the flattened poles are further from the source-containing sphere on which the 
spectrum is evaluated. Only at very high degree (about 720 for a terrestrial flattening) can this effect be 
overcome by the elevation of weaker, nearly horizontal domain dipoles on the equatorial bulge. 

The effect of scattering the domain dipoles throughout a thick ellipsoidal annulus amounts to 
multiplying the spectrum expected for a thin spherical shell by a thick ellipsoid factor TE n (for either type I 
or type II dipoles). The factor differs but slightly from the product of thickness and ellipticity factors. 

7.3 Comparison of Theoretical and Observational Crustal Spectra. 

When downwardly continued to the reference radius (6371 km), the observational magnetic spectrum 
R n derived from Magsat data by Cain et al. [1990] tend to increase with degree n for degrees 16 through 60. 
This conflicts with the notion of an exponential crustal spectrum that levels off at or below the Earth’s 
surface. Following Cain et al. [1989b], the increase in R n with n may be due to downward continuation of 
a noise component that gives a level spectrum at satellite altitude. The increase may also be due to, and is 
expected from, a magnetic field originating in randomly magnetized domains, each acting as a random 
dipole, scattered throughout Earth’s crust. When both the satellite altitude noise correction of Cain et al. 
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[1989b] and the spectrum expected from a shell of random dipoles are used to describe R n for degrees 16- 
50, the standard deviation of the residuals is minimal for an apparent source depth in the 20-55 km range. 

Although the covariance needed to weight R n was not available, higher degrees can be de- 
emphasized by considering satellite altitude R n (a + 420 km) or by trunction at degree 40. So doing tends to 
lower a source shell of random dipoles to seemingly implausible depths of 60 km or more (with or without 
the noise correction). It also reduces the leveling altitude of a plain exponential spectrum fitted to noise 
adjusted R n . We have yet to encounter a physical reason to suppose a level R„(r*) spectrum describes 
either a geophysical source, measurement error, or modeling error at radius r*. 

The results show that, in the absence of reliable uncertainty estimates, neither an elementary 
functional form for the crustal spectrum, a source shell depth, nor even a simple leveling depth can be 
reliably extracted from the field models considered. There is, however, some indication that a priori 
specification of a spectral form An a (r x /a) 2n ’ 2 with source depth a-r x in the 0-60 km range would return a 
power law index a: 0 < a < 3. Indeed, analysis of an equivalent source model of Magsat anomalies 
suggests a crustal refemce spectrum that varies approximately as 

7.4 Conclusion 

A suitable geologic null-hypothesis for the magnetic field orginating in Earth’s crust is that it results 
from an ensemble of uncorrelated domain dipoles, either of random orientation or random polarity, 
scattered throughout the lithosphere. The spatial magnetic power spectrum due to such a randomly 
magnetized crust differs from the exponential form describing geometric attentuation with radius above a 
sphere on which the spectrum is independent of spherical harmonic degree n. For harmonic degrees 16-60 
of interest rn satellite anomaly studies, the difference amounts to multiplying the exponential spectrum by a 
factor of n . Mathematically interesting corrections for non-zero uniform crustal thickness and ellipticity 
are of little or no practical concern in this range; they are important for harmonic degrees approaching the 
reciprocal relative thickness, or the reciprocal flattening, of the magnetized layer. 

Analysis of observational spectra from Magsat by Cain et al. [1990] on the reference sphere indicates 
that the spectrum expected from random dipoles on a shell offers more satisfactory perfoimance that a plain 
exponential spectrum. The former gives a lithospheric source depth for crustal fields instead of an 
ionospheric source altitude. This is also true when the reference spectrum is adjusted for a satellite altitude 
noise level. However, analysis of the noise adjusted observational spectrum at satellite altitude indicates 
that the random dipoles spectrum gives excessive source depths, while the plain exponential spectrum gives 
a broad range of leveling altitudes. The contrast suggests that the dominant part of the crustal magnetic 
spectrum is of an intermediate form which neither levels off, nor increases as rapidly as n 3 , on a 
lithospheric source shell. Such a form should be expected from correlated domain dipoles describing 
coherently magnetized geologic structures. 
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APPENDIX A 

At the top of the viscous sub-layer, regarded as the sphere of radius c, the well-known magnetic 
induction equation. 


3 t B = Vx(vxB) + (po)' 1 V 2 B , 

reduces to the diffusion equation because of the no-slip boundary condition on the fluid velocity v(c,0,<t>;t) 
= 0. Here p and a respectively denote core magnetic permeability and core electrical conductivity at the 
top of the sub-layer. The diffusion equation for the radial component of the core field is 

d t B r = (Moc)- 1 V 2 rB r = (pa)' 1 /d T \ , (A2) 

where /5 r B f / denotes the jump in the second radial derivative of B r across the core-mantle interface. This 
jump in the curvature of the field lines is due to the jump in lateral electric current density across the 
boundary between the conducting core and the resistive mande. Following Voorhies & Conrad [1996], the 
spherical harmonic expansion of (Al) is summarized via its coefficients 


dtgj^Kn+iX^c)”* 2 ] s (^ic) _1 J n m 

s (fto)- 1 (Kkm m ) 2 g n ™ [(n+l)(a/c) n+2 ] . (A 4) 

Equation (A3) defines the current density jump coefficient J n m ; (A4) defines an effective modal scale 
height [k nj m (t)] 1 which may be either real or imaginary. Similar relations can be written for cf,h m in 
terms of J^ m and k^ m (t). 

Summing the squares of these relations over orders m and dividing by (n+ 1) yields 


F nc(c) = (* 2 /pa) 2 (n+ 1 )(a/c) 2n+4 


1 t lk m m ) 4 (g n ?) 2 + 0^ m ) 4 (hnc > 2 ] 


(A5) 


= (po)' 2 (Ttk^) 4 R nc (c) (A6) 

which defines positive k^ 4 [k rn ] _1 is the effective scale height for harmonic degree n. Upward 
continuation of (A5) to any radius r > c gives 


F nc = (HO)' 2 (7tkm) 4 R nc . ( A7 ) 

The expectation value of (A7) yields (4c), with the understanding that the contribution to (F n ) from 
lltk m) l R n* = ( rek m^ 4lR nJ is of more immediate concern than are correlated fluctuations in (k t[| ) 4 and 
R n- 


APPENDIX B 

The integral in the ellipticity correction (35) is 



where 

a 2 * (y 2 + e 2 )/(l - e 2 ) = ^(2 - e 2 )/(l - e 2 ) 


(Bla) 
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(Bib) 


= 2y^ + V 4 . 

Because y 2 « 1, we approximate (1 + aV) 1 / 2 = 1 + y2 x 2 an(J 

Qn “ / () (l + 7 2 x 2 ) n dx = ZJy. (B2) 

The integrand is expanded in powers of ( y x) 2 « 1 


,1 00 (n+j-1)! 0 j 

Z„ =7/ Z - — — — - (-7V x 2 J dx , 

" 0 j=0 (n-l)!j! 

integrated term by term, and evaluated at the limits 

z _ £ (n+j-1)! C-^ 1 ) 
n " j=0 (n-l)!j! (2j+l) 

Equation (B4) is rewritten as the sum of indefinite integrals 
n (n+j-1)! . 

Z n = Z T / ('Y) 2j <*Y 
j =0 (n-l)'j' 

which is equal to the indefinite integral of the sum 
. n (n+j-1)! 0 . 

Z„ = J 2 — — (- 7 T 2 ) J dc 


j=0 (n*l)U! 

= J(1 + y 2 )' n dy . 


To obtain the recursion relation for Z n , note 


(B3) 

(B4) 

(B5) 

(B6) 

(B7) 


— fY(l + Y 2 )'"] = (1 +7 2 r n - +7 2 )' n ‘ 1 (B8) 

dy 

= -(2n - Ifr^l + T 2 )* 0 ' 1 + (1 + t 2 )'"' 1 ; (B9) 

therefore, integration of (B7) by parts and back substitution from (B9) yields 

z n * Y(1 + Y 2 )’ 11 + 1n\y*- (1 + 7 2 )'"' 1 dy (BIO) 

- ( 2nZ n+l * Y(l+Y 2 )’ n ]/(2n-l). (Bll) 

Equations (Bll) and (B2) give the recursion relations 

Qn+1 - [(2n-l)Q n +(l + y 2 )- n ]/(2n). (B12) 


Detailed calculations, inititalized with Qj = 1 - y 2 ^, show that (B12) underestimates Q n by about 4 parts 
per million. This is mainly due to ommission of a factor of about (1 + ytyl2). 
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